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ABSTRACT 

The cluster Praesepe (age ~ 650 Myr) is an ideal laboratory to study stellar evolution. 
Specifically, it allows us to trace the long-term decline of rotation and activity on the 
main-sequence. Here we present rotation periods measured for five stars in Praesepe 
with masses of 0.1 — 0.5 Mq - the first rotation periods for members of this cluster. 
Photometric periodicities were found from two extensive monitoring campaigns, and 
are confirmed by multiple independent test procedures. We attribute these variations 
to magnetic spots co-rotating with the objects, thus indicating the rotation period. 
The five periods, ranging from 5 to 84 h, show a clear positive correlation with object 
mass, a trend which has been reported previously in younger clusters. When comparing 
with data for F-K stars in the coeval Hyades, we find a dramatic drop in the periods 
at spectral type K8-M2 (corresponding to 0.4 — 0.6 Mq). A comparison with periods of 
VLM stars in younger clusters provides a constraint on the spin-down timescale: We 
find that the exponential rotational braking timescale is clearly longer than 200 Myr, 
most likely 400-800 Myr. These results arc not affected by the small sample size in 
the rotation periods in Praesepe. Both findings, the steep drop in the period-mass 
relation and the long spin-down timescale, indicate a substantial change in the angular 
momentum loss mechanism for very low mass (VLM) objects, possibly the breakdown 
of the solar-type (Skumanich) rotational braking. While the physical origin for this 
behaviour is unclear, wc argue that parts of it might be explained by the disappearance 
of the radiative core and the resulting breakdown of an interface-type dynamo in the 
VLM regime. Rotational studies in this mass range hold great potential to probe 
magnetic properties and interior structure of main-sequence stars. 

Key words: stars: low-mass, brown dwarfs, stars: rotation, stars: evolution, stars: 
activity 



1 INTRODUCTION 

Clusters are ideal environments to study the evolution of 
stars. Members of one particular cluster share age, distance, 
as well as metallicity, and these parameters are typically 
well-constrained. Thus, clusters provide 'snapshots' of the 
stellar population at a given age. Comparing the proper- 
ties of stars in clusters of different ages allows us then to 
track their evolution. Furthermore, clusters are often dense 
enough to allow multiplexed observations, i.e. wide-field 
imaging or multi-object spectroscopy, covering a large sam- 
ple of objects simultaneously. 

Two key parameters of stars, which strongly depend on 
age, are rotation and magnetic activity. The fundament for 
the understanding of the main-sequence evolution of rota- 
tion and activity is the empirically found Skumanich law: 

* E-mail: asllO@st-andrews.ac.uk 
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Looking at averages of large samples of objects, the rota- 
tional velocities of F to K stars evolve proportional to the 
inver se squareroot of age (e.g., ISkumanic hi 1 19721 : iBarnesI 
An analogous decline is seen for indicators of mag- 
netic activity, implying a correlation between rotation and 
activity, which has b een confirmed in mult iple studies (e.g., 
IStauffer et all Il997l ; [Terndrup et ai1l2000t ). This is usually 
interpreted with a twofold connection between stellar ro- 
tation and magnetic activity, providing a strong coupling 
between those two parameters: a) The interior dynamo of 
F-K stars, i.e. the underlying mechanism to generate mag- 
netic activity, is thought to be governed by rotation, b) The 
rotational braking on the main-sequence is due to angu- 
lar momentum losses th rough magnetic stellar winds (see 
ISchriiver fc Zwaadl200oh . 

To explain the presence of large-scale, stable magnetic 
fields in solar-type stars, it is usually assumed that the dy- 
namo in F-K stars is similar to the solar-type au> dynamo 
and operates in the transition layer between the convective 
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envelope and the radiative core. Therefore, it is expected 
that the magnetic field generation, and thus the magnetic 
and rotational properties, change when objects become fully 
conve ctive at masses < 0.3 — 0.4 jChabrier fe Baraffd 
1 19971 ). However, no clear change of rotation and activity in- 
dicato rs has yet been dete cted at or around this mass limit 
(e.g., iDelfosse et al. I ll998h . In addition, it is unclear what 
type of magnetic field is harboured by fully convective ob- 
jects, and if it depends on rotation. Investigating the ro- 
tational evolution and angular momentum loss of very low 
mass (VLM) objects on the main-sequence can contribute 
to clarify these issues. 

An ideal laboratory to analyse the rotation of VLM 
main-sequence stars is Prae sepe: With an age probably be- 
tween 0.6 and 0.9 Gyr (see iMermilliodl 1 198 J : I Adams et all 
2002), the cluster is significantly older than zero age main 
sequence (ZAMS) clusters like the Pleiades or a Per. At 
these ages, all stars have safely arrived on the main-sequence 
and the rotational regulation is thus dominated by magnetic 
wind braking, and not anymore affected by initial condi- 
tions, disk-braking, and contraction - the objects have es- 
sentially 'forgotten' their pre-main sequence history. Thus, 
it is possible to isolate the effect of wind braking on the 
rotational properties. 

Compared with the similarly old clusters Hyades and 
Coma, P raesepe is relatively co mpact, yet not too far away 
(180 pc, iRobichon et al.lll999T ). allowing us deep, efficient 
observations with wide-field facilities. Here we report on 
photometric monitoring for 18 members of Praesepe, which 
provided the first five rotation periods for members of this 
cluster - all for objects with estimated masses below 0.5 Mq. 
Combine d with previously publis he d peri ods for younger 
objects (|Scholz fc Eisloffell l2004al lbl l2005h . this allow us 
to probe the long-term evolution of rotation in the VLM 
regimeQ 



2 PHOTOMETRIC MONITORING 
2.1 Observations 

The analysis in this paper is based on two photometric time 
series of VLM stars in Praesepe. The observing strategy in 
both runs was similar: To mimimize internal inconsistencies 
in the data and thus systematic uncertainties, we 'stared' at 
a particular field in Praesepe. We did not 'dither' around 
a central position and made an effort to position the field 
of view at the same coordinates in each observing night. 
Therefore each object was at roughly the same pixel posi- 
tion on the detector throughout the run. The exposure time 
for the single frames was 600 sec in both runs. Since we were 
interested in very low mass objects, which have a spectral 
energy distribution peaking in the near-infrared, all obser- 
vations were carried out in the I-band filter. 

The target fields were selected to maximise the num- 
ber of known Praesepe members in the field of view. We 



1 Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck Institut fur Astronomic and the Institute de As- 
trofisica dc Andalucia (CSIC), and at the Thiiringer Landesstern- 
warte Tautenburg (Germany) 



Table 1. Logfilc of the time series observations: Run, date of 
observations, no. of images, typical seeing 



Run 


Date 


No. 


Seeing 


TLS 


16/01/2001 


12 


2'.'0 




17/01/2001 


12 


2'.'5 




18/01/2001 


6 


2"0 




14/02/2001 


32 


2"0 




15/02/2001 


38 


1'.'8 




16/02/2001 


21 


l'.'S 




18/02/2001 


1 


2" 5 


LAICA 


23/01/2003 


13 


3'.'8 




24/01/2003 


19 


2'.'5 




25/01/2003 


11 


2 / /0 




27/01/2003 


22 


1'.'5 




28/01/2003 


43 


1'.'3 



compiled a n initial catal o gue of VL M members f r om th e 
surveys of iHamblv et all (|l995al lbh: iPinfield et ail l|l997h : 
iMagazzu et al.l ( 1998h . at the time of the observations (2001, 
2003) essentially the complete census of the VLM population 
in this cluster. For all these objects, the main criterium for 
cluster membership is multi-colour photometry, confirmed 
in many cases by proper motion s and /or spectroscopy . Ac- 
cording to a follow-up study bv lHodgkin et ai] (|l999h , the 
contamination rate is ~ 10% in the Hambly et al. sample 
and ~ 50% among the (fainter) Pinfield et al. objects. 

The first time series was obtained using the 2-m 
Schmidt telescope at the Thiiringer Landessternwarte Taut- 
enburg (TLS, Germany), equipped with a 2048 x 2048 SiTe 
CCD camera. The wide-field camera provides a spatial reso- 
lution of l'.'2/pixel, resulting in an unvignetted field of view 
of about 0.36sqdeg. The observations at the TLS were car- 
ried out in two sessions in January and February 2001, with 
a gap of almost four weeks. Thus, although we cover more 
than a month in total, our sensitivity to long periods (1-4 
weeks) will be limited. The observations were partly affected 
by cirrus and mediocre seeing conditions. In total, the run 
provided a time series of 125 images of the same field. The 
target field for the TLS run was centred at a = 8 h 36 m 53 s , 
5 = +19°48'24" (J2000.0). This field contains nine mem- 
bers from the Hambly et al. catalogue (H85, 91, 102, 106, 
115, 126, 140, 158, 181 in their nomenclature) and two from 
the Pinfield survey (PI and P2 in their nomenclature), plus 
some more higher-mass members which are saturated in our 
images. 

A second campaign was carried out using the wide-field 
imager LAICA at the 3.5-m telescope at Calar Alto Obser- 
vatory from 23-28 Jan 2003 (in service mode). LAICA is 
a 2 x 2 mosaic of 4 x 4K CCDs, which we used to 'stare' 
at our target field in Praesepe. Thus, we did not observe a 
continuous field, but four 'patches' of 15' x 15', separated 
by gaps of similar size, in total 0.25 sqdeg at a resolution of 
O'.'22/pixel. The total length of the campaign was 10 nights 
(scheduled as 10 half nights), but most of the first half of the 
run was not useful due to excessively high seeing and cloud 
coverage. In the last six nights, however, we were able to 
obtain in total 108 images with mostly good quality. Thus, 
this time series will be highly sensitive to periods up to a 
few days. The LAICA run was focused on a target field at 
a = 8 /l 40 m 56.5", S = +19°32'30" (J2000.0), not overlapping 
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with the TLS field. Since the LAICA images are consider- 
ably deeper than the TLS data, we aimed to cover some of 
the faintest member candidates in the cluster. The field thus 
contains six objects from the Pinfield survey (P14, 15, 16, 
17, 19, 20) plus one additional unsaturated star from the 
Hambly list (H218). 

2.2 From images to lightcurves 

Image reduction, photometry, and relative calibration for 
both runs followed the routines established in the fram e- 
work of our previous campaigns (|Scholz fc E isloffcl 200lallbl. 
2005). In the following, we give a brief account of the 
basic principles of our method to derive lightcurves from 
time series images. For details, we refer to the more elabo- 
rated discussion in the af ore mentioned papers, particularly 
IScholz fc Eisloffell | 2004al ) . In case of the LAICA campaign, 
all reduction and calibration steps were carried out sepa- 
rately on the four individual chips. 

The image reduction started with bias subtraction and 
flatfield correction (based on high signal-to- noise domeflats). 
To remove the interference ('fringe') pattern produced by 
night sky line emission, which is typical for deep I-band 
exposures, we constructed a fringe mask (an image of the 
fringes, without objects) by averaging a number of dark sky 
exposures. This fringe mask was then properly scaled and 
subtracted from the individual time series images. Scaling is 
necessary because the amplitude of the fringes is a function 
of airmass and observing conditions. After this process, the 
residuals of the fringes in the final images do not exceed 1% 
of the sky background. 

Using the Source Extractor, we created an object cat- 
alogue for one selected image in each time series ('master 
image'). Customised IRAF routines were used to a) deter- 
mine pixel offsets between the master frame and all individ- 
ual time series images, b) creating object catalogues for all 
time series images, c) carry out aperture photometry for all 
objects. For the TLS run, PSF fitting based on daophot was 
chosen as the optimum method for photometry. We selected 
a sample of PSF reference stars and used them to model 
the PSF in all time series images. Based on the values from 
the aperture photometry, new magnitudes were derived by 
adapting the model PSF to all objects in the field. 

We refrained from using PSF fitting techniques for the 
LAICA run: Due to the strongly variable seeing conditions, 
it was not possible to find a set of PSF reference stars usable 
for a complete time series. Changing the reference PSF, how- 
ever, may introduce systematic uncertainties. On the other 
hand, LAICA has an excellent pixel resolution, diminish- 
ing the benefits of PSF fitting photometry in moderately 
crowded fields. Therefore, we rely on the more robust aper- 
ture technique. 

The instrumental magnitudes are still affected by the 
combined effects of varying airmass, varying atmospheric 
extinction, and variable zeropoint. This is corrected by sub- 
tracting a 'master lightcurve' from all time series, which is 
an average lightcurve calculated from a set of non-variable 
reference stars. The selection of these reference stars is the 
critical part in the relative calibration of the photometry. In 
a first step, we chose only objects with valid photometry and 
photometric error below 0.1 mag in all images. From this ini- 
tial set (typically a few hundred objects) , we excluded poten- 



Table 2. Objects with si gnificant periodic vari a bility in Prae- 
sepe: Object ID (following THamblv et ail jl995bT) ; [Pinfield et all 
(1997)), mass estimate (see Sect. l4.il . spectral type estimate (see 
Sect l4.ll l. period, period uncertainty, period amplitude 



Object M (M Q ) Sp.T. P (h) AP (h) A (mag) 



TLS: 












P2 


0.12 


M6.5 


5.65 


0.01 


0.04 


H91 


0.29 


M3 


42.2 


0.20 


0.08 


H115 


0.40 


M2 


83.6 


1.02 


0.04 


LAICA: 












P20 


0.11 


M7 


12.2 


0.20 


0.07 


H218 


0.19 


Ml 


16.2 


0.51 


0.04 



tially variable objects by testing the rms in their lightcurve 
against the rms of all other preliminary reference stars. In 
several steps, the sample of reference stars is cleared from 
contamina ting variable objects. For details of the process, 
we refer to IScholz fc Eisloffell (1200481 ) . After subtracting the 
resulting master lightcurve from all time series, we obtained 
lightcurves in relative magnitudes for all objects in our fields. 
This database is the fundament for the variability analysis 
described in Sect. [3] 

We give two estimates for the accuracy of the relative 
calibration: a) The average rms of the final lightcurves of the 
reference stars (found to be non-variable) is 0.01 mag for the 
TLS and 0.02 mag for the LAICA run. b) The minimum rms 
measured from the lightcurves of bright stars is 0.009 mag in 
the TLS run and 0.01 mag for the LAICA run. Both methods 
provide an estimate of the intrinsic noise in the lightcurves. 
We conclude that for bright, but unsaturated objects we 
achieve a photometric accuracy of about 1% in both runs. 



3 PERIOD SEARCH 

The main goal of the photometric monitoring was to de- 
tect periodic changes in the flux of our targets, which can 
be interpreted as indication of co-rotating spots - thus to 
measure rotation periods. The lightcurves of the Praesepe 
members in our fields were scrutinised with a rigorous pe- 
riod search procedure, which has been implemented and 
extensively tested in the fr amew ork of previous campaigns 
jScholz fc Eisloffell l2004a1 lbl, 120051 ). 

Period search in astronomic time series is non-trivial, 
mainly because of the characteristic gaps caused be day- 
time, weather changes, and unavailibility of the telescope. 
As a result, the distribution of datapoints is usually strongly 
clumped, hampering time series analysis. This is further 
complicated in the case of very low mass objects, where 
signal-to-noise (defined as ratio between amplitude of the 
periodic variation and noise level in the lightcurve) is typ- 
ically less than five, so that the period is not always obvi- 
ous from a by-eye examination. We try to mitigate these 
problems by a) maximising the number of datapoints in the 
lightcurve (to obtain as much information as possible about 
the light changes) , and b) by probing for periodicities using 
a variety of independent and robust tests. In the following, 
we briefly outline the criteria which we require to be fulfilled 
to accept a period: 

(i) The Scargle periodogram (jScarglel Il982h shows 
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a peak with a preliminarily determined false alarm 
probability (FAP) belo w 1% (FAP calculated following 
iHorne fc Baliunasl (|l986l )). 

(ii) Nearby reference stars do not show the same peak in 
the periodogram. 

(iii) The scatter in the lightcurve is significantly (based 
on a F-test) reduced after subtraction of the sine-wave ap- 
proximated period. 

(iv) The lightcurves of reference stars phased to the de- 
tected period do not show an y sign of periodicity. 

(v) The CLEAN algorithm dRoberts et al.ll 19871 ) confirms 
the periodogram peak, i.e. it is not an artefact from the 
window function. 

(vi) The empirical FAP determined using a bootstrap ap- 
proach (see below) is below 1%. 

(vii) The phaseplot shows the period clearly (but see be- 
low). 



P2 



P = 5.65h 



H91 



P = 42.3h 



Applying these criteria, we identified five objects with 
significant photometric period in the lightcurve (see Table 
[2]). The phaseplots for these five objects are shown in Fig. 
[T] Please see below for individual notes on the lightcurves. 
Final values for the FAP have been det ermined using the 
bootstr ap approach explained in detail in lScholz fc Eisloffell 
(|2004al V The basic idea is to produce 10000 randomised 
lightcurves per object by retaining the observing times and 
shuffling the data-values only. The resulting lightcurves will 
have the same sampling and noise level as the original time 
series, but certainly no periodicity. For all 10000 random 
lightcurves we calculated the Scargle periodogram and de- 
termined the amplitude of the highest peak. The fraction 
of lightcurves, for which this number exceeds the peak am- 
plitude for the measured periodicity of the given object is 
a robust and reliable estimate for the FAP for this period. 
For all five detected periods, the peak in the periodogram 
is clearly higher than any peak from the 10000 randomised 
lightcurves; thus the FAP is below 0.01%. 

We tested the range of periods for which we are sensi- 
tive by adding sinewaves with typical amplitudes of 0.05 mag 
and varying period to the lightcurves of non-variable refer- 
ence stars. By applying the period search procedure used 
for our target stars to these 'artificial' periodic objects, we 
find that our sampling allows us to reliably detect periods 
between ~ 1 h and ~ 4 days (in both runs) . This is in line 
with the expectations: Both monitoring campaigns provided 
dense sampling over at least three consecutive nights, with 
only spotty coverage for longer timescales. Hence, it is not 
surprising that all our five detected periods are in the range 
between 1-4 d. However, it has to be emphasised that the 
period search is not sensitive for periods around multiples 
of 1 day, due to the daytime gaps, thus the period sample is 
probably not complete. On the other hand, both time series 
have some potential to detect periods longer than 4 days (up 
to periods of a month for the TLS run). Still, given the fact 
that our sensitivity for these long periods is limited, we do 
not put too much emphasis on the fact that we did not find 
any such period. 

Our period search still includes one subjective cri- 
terium: the by-eye examination of the phased lightcurves. 
We require that the period is visible in this plot, i.e. 
that it shows a clear maximum and minimum. However, 
this criterium has to be used with caution. As we have 



H115 



P = 83.6h 



0.5 

P20 



1.5 2 

P=12.2h 



— - ^ 



H210 



P=16.2h 



Figure 1. Phase plots for the five objects with significant periodic 
variations. P2, H91, H115 have been observed in the TLS run, P20 
and H218 in the LAICA run . © 2002 RAS, MNRAS 000,[JJflO] 
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shown in lScholz fc Eisloffel (2005), periodogram techniques 
are able to pick up periods even at signal-to-noise levels 
too low to produce a visually obvious periodicity in the 
lightcurve. Thus, a not entirely convincing phaseplot (see 
Fig. [TJ does not necessarily indicate a false detection. This 
has to be taken into account when discussing individual ob- 
jects. Please note that the objects P2, H115, and P20 have 
been con firmed as cluster me mbers by near-infrared pho- 
tometry (jHodririn et al.lll999l ). Based on their position in 
the colour-magnitude diagram, P2 and H115 have been clas- 
sified as possible binaries. 

Notes on individual objects: 

P2: The object is one of the faintest, lowest mass mem- 
bers of Praesepe discovered so far. We estimate roughly 
0.1 Mq; if the object really is a binary, this corresponds to 
the total system mass. Additionally, the amplitude of the pe- 
riodicity if fairly low, resulting in low signal-to-noise ratio. 
This is a typical case where the effects of the low signal- 
to-noise are offset by a large number of datapoints (in this 
case 123). Scargle and CLEAN periodogram show a convinc- 
ing peak, which translates into a highly significant period of 
5.65 h. 

H91: Although the phaseplot for the period of 42.2 h ex- 
hibits some gaps, the datapoints do show a clear flux mod- 
ulation, i.e. maximum and minimum are clearly visible. 

H115: Due to the relatively long period of 83.6 h, the 
coverage is not continuous. The two strong clumps at phases 
of 0.3 and 0.85 represent the datapoints from Feb 14 and 16, 
respectively, which are clearly offset by roughly 0.04 mag, 
and thus give direct evidence for a photometric variation on 
timescales of days. 

P20: The phaseplot shows a clear periodicity, but the 
lightcurve appears 'box'-shaped. This raises the question 
whether we are dealing with an eclipsing binary rather 
than a rotational modulation. However, when examining 
the lightcurve, there are clear smooth and gradual trends 
in some nights rather than rapid brightness changes, which 
is difficult to explain in an eclipse scenario. Moreover, the 
fact that the suspected eclipse covers almost half of the 
phase space would require a roughly equal-radius binary sys- 
tem, which is unlikely given the very small eclipse depth. 
Moreover, the period would imply an orbital separation 
of fractions of 1 AU, a nd such systems seem to be rare 
jFischer fc Marcvlll992T ). Therefore, we attribute the 'boxy' 
shape of the lightcurve to an unusual spot configuration. 
Still, the possibility of an eclipsing system cannot completely 
be disregarded. 

H218: The largest peak in the Scargle periodogram cor- 
responds to a period of 33.6 h, but disappears after apply- 
ing the CLEAN routine and is thus probably an artefact. 
CLEAN, however, identifies a peak at the doubled frequency 
(corresponding to a period of 16.2 h) as real. Since the same 
peak is also highly significant in the Scargle periodogram, we 
adopt a period of 16.2 h for this object. The period remains 
significant when we exclude the five outlying datapoints at 
phases of 0.15 and 0.7. 



4 INTERPRETATION 

The best explanation for the observed periodic variability 
is the presence of spots on the surface of the stars, co- 



rotating with the objects and thus modulating the flux. With 
an idealised spot distribution (one spot, regularly shaped), 
we would expect sine-shaped periodicities, particularly at 
low signal-to-noise ratio. In contrast, other sources of peri- 
odic variations (eclipsing binaries, planetary transits) pro- 
duce regular dips in the lightcurves. Four out of five periods 
clearly resemble a sinecurve. The sole exception, P20, could 
be an eclipsing binary, but we argue that its box-shaped 
lightcurve is more likely caused by a particular distribution 
of surface features. 

The underlying reason of the surface features is most 
likely magnetic activity: Our targets have effective tem- 
peratures larger than 3000 K (spectral types around mid 
M), rendering the possibility of condensated dust clouds, 
as often discussed for the much cooler L dwarfs, improb- 
able. M dwarfs, on the other hand, are known to har- 
bour substantial magnetic fields, as evidenced by X-ray and 
Ha activity. Specifically, for three of our periodic objects 
(H91, H115, H218) chromospheric activity evidenced by Ha 
emission with equi valent widths between 3 a n d 23 A has 
been reported by Barrado v Navascues et ail l|l998h and 



iKafka fc HonevcuttJ <|2006f ). It seems reasonable to assume 
that the remaining two objects (without Ha observation so 
far) share similar chromospheric properties. 

Thus, for the following discussion, we assume that the 
most likely origin for the observed periodic variations is a 
rotational modulation due to magnetically induced surface 
spots. Hence, the periods correspond to the rotation peri- 
ods of our targets. These are the first rotation periods mea- 
sured for objects in this mass range in a cluster significantly 
older than the ZAMS stage (which is reached at roughly 0.2- 
0.3 Gyr for very low mass stars, see llrwin et alj (|2007 )L In 
the following subsections, we will use this small period sam- 
ple in combination with literature data to probe rotation vs. 
mass and rotation vs. age for very low mass objects on the 
main-sequence. 

To assure ourselves that the five datapoints give a 
realistic first glance of the VLM period distribution in 
Praesepe, we carried out two plausibility checks. First, we 
compare with v sini values for the Hyades, compiled by 
iTerndrup et~aj] (|200Ch . The Hyades are roughly coeval with 
Praesepe, and since the v sin i sample size in the very low 
mass range is about 20, these values can provide a use- 
ful consistency check. In the mass range 0.2-0.4 Mq, the 
lower limit in rotational velocities is defined by the detec- 
tion limit of 5kms _1 , which translates to an upper limit 
in P sin i of 30-90 h , depe nding on mass (with radii from 
IChabrier fc Baraffd i| 19971 )1. The upper limit in v sin i is 
around 15kms _1 , corresponding to a lower limit for Psini 
between 10 and 30 h. Considering that sin i has a fairly tight 
distribution around 0.7-0.8 for random orientations of rota- 
tional axes, these numbers are completely consistent with 
our period limits. 

A second check can be made based on the period 
ranges for clusters somewhat younger than Praesepe. In 
the Pleiades, NGC2516, and M34, all clusters with ages be- 
tween 100 and 200 Myr, a substantial sample of VLM periods 
has been establis hed, in total ~ 3 0, the dominant major- 
ity in NGC2516 (llrwin et al.ll2006l . 120071 ; IScholz fc Eisloffell 
l2004bh . These periods range from 3 h to about 3 days. Given 
the fact that there is no known mechanism to accelerate ro- 
tation on the main-sequence, instead, we expect rotational 
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braking and thus longer periods in older clusters, this again 
seems consistent with the period range in Praesepe. A more 
detailed analysis of the rotational evolution based on a com- 
parison of Praesepe with younger clusters will be given in 
Sect. 14.21 Here we conclude that the small sample in Prae- 
sepe is likely to represent a typical range of VLM rotation 
periods in this cluster. 



4.1 Rotation vs. mass 

It is well-estab lished that rotatio n is a function of ob- 
ject mass (see iHerbst et all 120071 . for a review). Specifi- 
cally, there is accumulating evidence that in the very low 
mass regime the average period drops steadily with decreas- 
ing mass. This positive period-mass correlati on at masses 
< 0.3 - 0.4 M p has been found in the ONC l|Herbst et all 
l200ll). eOri ilScholz fc Eisloffell 120051 ), and the Pleiades 
1 Scholz fc Eisloffell 12004^ 7 clusters with ages between 1 
and 125 Myr Sim i larly in the clusters NGC2516 and M34 
l|lrwin et al.1 120061 . 120071 ). which roughly mark the ZAMS 
for very low mass objects at ages of 150 and 200 Myr, there 
seems to be a general decline of the upper enve lope of the pe- 
riods with decreasing mass (jlrwin et al.ll2007l . see their Fig. 
17). It has been pointed out that this trend is consistent 
with constant an gular momentum for all object masses (e.g. 
IHerbst et al.l [2001T ). and thus might just reflect the drop in 
stellar radius. Since it is already established at very young 
ages, it may be related to the initial distribution of angular 
momentum. 

Here we probe if our small period sample in Praesepe 
- older than all other clusters with rotation periods in this 
mass range - allows us to see a similar trend. We derived 
masses for the five objects in Table [2] by comparing the 
available photometry in t he I-band from the literature with 
evolutionary tracks from iBaraffe et al.1 (|l998l ) for an age of 
0.75 Gyr. According to these estimates, all five objects are 
in the very low mass regime with masses ^ 0.4 Mq. Due 
to age uncertainties, photometric band inconsistencies, and 
model shortcomings, the uncertainties in the derived masses 
are considerable (probably on the order of 50%). However, 
most of this uncertainty is systematic, thus we expect that in 
a relative sense our masses are realistic. We caution against 
comparing these values with masses derived using a different 
approach. Using the same model isochrone, we determined 
effective temperatu res for our five target s. Comparing with 
the T e ff scale from iLuhman et al.1 (2003) gives an estimate 
for the spectral types (see Table [5J ■ The uncertainties in 
these 'photometric' spectral types are probably ±1 — 2 sub- 
classes. 

In Fig. [2] upper panel, we plot periods vs. masses for 
the five objects in Table [2] It is immediately obvious that 
the periods appear to increase with mass in a roughly linear 
way. A linear least-square fit (shown as dotted line) gives: 
P = (240 ± 40) (M/Mq) - (21 ± 10) h. The correlation is 
significant with a false alarm probability of 1.0%. The slope 
in the relationship is clearly steeper t han in the Pleiades 
(105 ± 61 (M/M ), IScholz fc Eisloffell (|2004bl l). maybe in- 
dicating that rotational braking on the main-sequence is a 
function of object mass. Based on only five datapoints, how- 
ever, this finding is of somewhat limited value. Clearly, more 
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Figure 2. Upper panel: Rotation period vs. mass for VLM ob- 
jects in Praesepe. The solid line is a linear fit; the dashed line 
shows the P(M) fit for ou r VLM period sample in the Pleiades 
l|Scholz fc E isloffel 2004b). Lower panel: Rotation period vs. 
spectral type (as indicator for mass): Crosses are the five VLM 
periods in Praesep e from this work, tr iang les are periods in th e 
coeval Hyades from lRadick et all (l987) and Prosser et al.Nl995t) . 
Spectral type is parameterised as follows: GO - 10, K0 - 20, M0 
- 30. 



datapoints are needed to solidify the main-sequence P-M 
correlation in the VLM regime0 

It is more instructive to look on the VLM periods in 
comparison with rotation measurements obtained for more 
massive stars. To our knowledge, such data is not available 
for Praesepe, but for the Hyades, with 600 Myr a ro ughly 
coeval cluster (|Mermilliodlll98ll : iPerrvman et al.lll99Sh . For 
this comparison, we prefer to use spectral types instead of 
masses, in order not to be biased by model inconsistencies 
in the mass estimates, which are difficult to avoid when 
covering a mass range of more than one order of magni- 
tude. Since we use only coeval objects, the spectral type 
is a valid indicator for stellar mass. We collected a sam- 
ple of 25 periods for Hyades members from the photomet- 



2 We note that for the three objects with published Ha equiva- 
lent widths the rotation periods seem to decline with increasing 
activity, so it may be that the observed period-mass trend is in 
fact a reflection of activity levels increasing with later spectral 
types. 
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ric monitoring campaign s published bv lRadick etU] (|l987t l 
and iProsser et all (J1995) . Spectral types for these ob jects 
have originally b een published by Ivan Buerenl (|l952h and 
van Altenal(ll969l ). Combined with our sample, the spectral 
range from late F to late M is covered, roughly correspond- 
ing to a mass range from 0.1 to 2 Mq. 

The period/spectral type relation is shown in Fig. [2] 
lower panel. It clearly demonstrates that for F-K stars the 
periods increase to wards later spectral types. According to 
iRadick et ail(|l987f ). these periods can be explained in terms 
of the correlation between magnetic activity and the inverse 
Rossby number Ro, the ratio between rotation period an d 
convective turnover timescale re (e.g. iNoves et al.lll984l ). 
This relation basically implies that the magnetic field am- 
plification mainly depends on convection properties and ro- 
tation - supporting the idea of an au dynamo operating in 
F-K stars. We note that the Hyades periods for F-K stars 
represent what iBarnes! (|2003h called the Fsequence of rota- 
tional evolution, i.e. objects with interface dynamo. 

As it is apparent from Fig. [2] this sequence breaks down 
roughly at spectral types K8-M2, corresponding to a mass of 
about 0.4-0.6 M (or a B - V colour > 1.5). Without clear 
transition regime, the periods drop by almost one order of 
magnitude between K8 and M2. Even if the period distribu- 
tion at very low masses is not complete (but see above), the 
simple fact that the fast rotators do exist among M dwarfs 
and not in the G-K spectral range indicates a fundamental 
change in stellar rotation at the given mass limit. 

The underlying physical reason for this breakdown is 
not clear. Obviously, it indicates a dramatic change in the 
rotational braking at a certain mass, which may be related 
to a change in wind properties, surface magnetic field, or 
dynamo action. It should be noted that the break in rota- 
tion period in Praesepe occurs roughly at the same spec- 
tral type as the onset of chromospheric activity measured as 
Ha em ission: As recently reported by Kafka fc HonevcuttJ 
(2006), Praesepe objects earlier than Ml are rarely chromo- 
spherically active, in contrast to later type objects. Since 
enhanced activity would normally be associated with the 
potential of efficient angular momentum loss, this poses the 
question why the fast rotating M-dwarfs in Praesepe are the 
most active objects. 

In the rotational evolution scheme bv lBarned (|2003h . the 
fast rotating M-dwarfs in Praesepe would represent the C- 
sequence (quote: 'possess only a convective field', 'no large 
scale dynamos'). In this view, the mass limit at which I- 
sequence and C-sequence bifurcate should be a function of 
age and shift to later spectral types as objects get older. 
There is indeed some evidence for this: In v sin i data for field 
M dwarfs, most of them certainly older than Praesepe, the 
fraction of fast rotators incr eases from zero at M3 to 100% 
at M6 l|Delfosse et al.lll998l ). clearly at lower masses than 
the transition seen in Praesepe/Hyades. Still, this does not 
ex plain why the m agnetic field and the dynamo, as suggested 
bv lBarned ll2003h. should change abruptly at a given mass. 

Is there a link between the rotational properties and 
the interior structure of the stars in the VLM range? Ob- 
jects with 0.5 Mq are still thought to har bour a substantial 
radiative core (|Chabrier fc Baraffd H997T I. Going to lower 
masses, however, the bottom of the convective envelope 
drops quickly until the objects are fully-convective. There 
is strong reason to believe that dynamo action changes in 



a fundamental way with the disappearance of the radiative 
core. The aui dynamo ('interface dynamo') thought to oper- 
ate in solar-type stars requires a shear layer at the bottom of 
the convection zone, which is not supposed to exist in fully 
convective objects. Thus, it can be expected that the F-K 
type activity vs. Ro relation breaks down as soon as the ra- 
diative core disappears. At the same time, the objects may 
lose the ability of efficient rotational braking, possibly due 
to a change from large-scale to small-scale fields, resulting 
in fast rotators. 

However, structural models of low-mass stars essentially 
agree on the fact that this transition is supposed to occur at 
masses between 0.3 and .4 M fl for M/H = -2 .0 (e.g. 
IChabrier fc Baraffd 1 19971: iMontalban et all feood ) or even 
below l|Mullan fc MacDonaldl 1200 ll ). i.e. at lower masses 
than the observed drop in rotation periods apparent in 
Fig. [2] This limit depends on metallicity in a complex 
way and is expected to i ncrease somewhat for M/H > 
llChabrier fc BaraffdH997h . Since both Hyades and Praesepe 
have higher than solar metallicities (M/H of 0.17 and 0.14 
respectively), this might account for some of the the differ- 
ence in object mass between observed and expected transi- 
tion to fast rotators. If the break in the mass-period relation 
is entirely to explain with the change to fully convective ob- 
jects, however, remains questionable. 

The appearance of Fig. [2] and its possible potential 
to probe magnetic and structural properties is intriguing 
enough to justify future research to explain its origin. The 
following approaches appear to be particularly useful: a) 
modeling the transition to fully convective objects for a 
range of metallicities, b) complementary observations in the 
K8-M2 spectral range to check for changes in magnetic prop- 
erties (e.g., Doppler imaging, magnetic field measurements, 
multi- wavelength monitoring), c) rotational studies for stars 
with masses between 0.6 and 0.3 Mq, which are so far only 
represented by 4 datapoints. 



4.2 Rotation vs. age 

The rotational evolution of stars is believed to be controled 
by mainly three mechanisms: 

a) disk braking, i.e. angular momentum losses due to cou- 
pling between star and circumstellar di sk, for example as s o- 
called 'disk-locking' (see the review bv lHerbst et~alll2007l ). 

b) pre-main sequence contraction resulting in a spin-up, 

c) magnetic winds driven by dynamo actions, carrying away 
mass and angular momentum. 

In Fig. 13] we compare our rotation periods for Praesepe 
with period data for VLM objects in younger clusters: eleven 
rotation periods for stars in the Pleiades (ag e 125 Myr), 
l|Scholz fc Eisloffelll2004bl : lTerndrup et al.lll999T ). about 250 
periods for objects NGC2516 ( age 150 Myr) show n as me- 
dian, 10%, and 90% percentil es (llrwin et alj|2007l ). nine pe- 
riods in M34 (age 200 Myr, llrwin et al.l ~([2006). Together 
with the five Praesepe periods, this is to our knowledge the 
total available period data for main-sequence VLM stars. 
In all four clusters, the sampling of the photometric mon- 
itoring was sensitive to periods ranging from 30min to 
at least four days. In the Pleiades and M34, the detec- 
tion range extends to 12 and ~ 10 days, respectively. To 
illustrate the differences in the mass ranges, objects with 
0.2 < M < 0.4 Mq are shown as plus signs and lower mass 
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objects as crosses. Note that the majority (~ 90%) of the 
periods in NGC2516 has been obtained for stars in the upper 
mass bin . For Praesepe, we use here an age of 650 Myr, as 
given by iMermilliodl (|l98ll ); shifting the datapoints within 
600-900 Myr (the likely age range for the cluster) does not 
affect the following results. 

We calculated simple rotational evolutionary tracks tak- 
ing into account spin-up due to contraction and angular 
momentum losses by stellar winds. Since typical disk life- 
times are < 10 Myr, the influence of disk braking is irrele- 
vant here. The spin-up was estimat ed using radii from the 
models of lChabrier fc Baraffd l|l997h . Since even the lowest 
mass objects have reached their final radii to within 10% af- 
ter 200 Myr, the effect of contraction of the rotation period 
is only visible between 100 and 200 Myr. After that, the ro- 
tational evolution is almost entirely determined by braking 
due to stellar winds - that's why main-sequence clusters are 
the ideal environment to test wind braking. For solar-type 
stars, wind brakin g results in th e so-called Skuma n ich law 
l|Skumanichlll972l ): P oc y/i. In IScholz fc Eisloffell <|2004bh 
we demonstrated that this braking law is not applicable in 
the pre-main sequence evolution of VLM objects, because 
it predicts periods > 4 days at ~ 100 Myr, while the upper 
VLM period limit at this age, reliably determined from the 
available period and v sin i data, is only ~ 2 days. Thus, a 
more moderate braking law has to be used, as it is expected 
for objects w ith saturated dynamo (|Terndrup et al.l [20001 ; 
Barnes 2003) or for objects with a strong concen tration of 
magnetic flux near the pole (jSolanki et al.l [19971 ). We use 
here an exponential form for the braking law P oc exp(t), 
as it is expected for a saturated dynamo. Given the lack 
of understanding for the underlying physics for this type of 
rotational braking, this should be treated as an ad-hoc solu- 
tion to provide moderate braking, rather than an accurate 
physical model. 

Taken together, the rotational evolution can be ex- 
pressed as: (Ph: initial radius, Rj: final radius, Pi: initial 
period) 



P f = a x (Ri/Rf) 2 x P, 



(1) 



with a = 1 for zero wind braking and a — exp(t/r) for 
exponential braking, where r is the braking (or spin-down) 
timescale. In Fig. [3] we plot both cases, no braking with 
dotted lines, exponential braking with dashed lines. The ro- 
tational evolution is shown for 0.4 Mq, starting at Pi = 40 h 
(upper limit in the Pleiades) and 0.1 Mq, starting Pi = 3h 
(lower limit in the Pleiades). The different values for Pi were 
chosen to reflect the P-M relationship discussed in Sect. 14. ll 
As can be seen from the tracks, changing the mass does not 
significantly affect model tracks. 

In the discussion of this plot, it is important not to be 
misguided by the small number of periods. In particular, it 
is essential to keep in mind that the Praesepe period range 
is probably incomplete (in contrast to the younger clusters). 
Thus, rather than reproducing upper and lower period lim- 
its, the rotational evolutionary tracks should simply be able 
to explain the existence of the periods measured in Praesepe. 

As can be seen in Fig. [3] tracks without any braking on 
the main sequence (dotted lines) have problems reproduc- 
ing the longest periods measured so far in NGC2516, M34, 
and Praesepe. Thus, some rotational braking on the main- 
sequence is likely occuring among VLM objects, in agree- 
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Figure 3. Rotation perio ds for 0.1 — 0.4 mass objects in 
Pleiades, N GC2516, M34 l lScholz fc Eisloffell l2004bl ; llrwin etafl 
120061 . 120071 ) , and Praesepe (this paper) . We chose different sym- 
bols for two mass bins: 0.2 - 0.4 Mq - plus (+), 0.1 - 0.2 M 
— cross (x). For NGC2516, we plot the median as octagon and 
the 10%- and 90%-percentiles as horizontal lines. Note that the 
upper period limits in Pleiades, M34, and NGC2516 are reliably 
determined, while the period range in Praesepe may be incom- 
plete. The dotted lines show the period evolution as expected for 
zero wind braking. Dashed lines include exponential wind braking 
with r = 600 Myr. 



ment with our findings in IScholz fc Eisloffell (|2004bh . From 
the tracks for exponential braking, we can rule out that the 
braking timescale r is shorter than 200 Myr: Such low values 
for r would imply that an object with 3 h rotation period 
at 100 Myr (lower limit in the Pleiades) has a period longer 
than 100 h at the age of Praesepe - thus not allowing for 
any of our five measured periods. The plausible range for 
the braking timescales is between 400 Myr and 800 Gyr (in 
Fig. [3]we show the tracks for r = 600 Myr). 

These results are mostly consistent with earlier find- 
ings for the spin-down t imescale in the VLM regime. In 
IScholz fc Eisloffell (|2004bl) we conclude that the rotational 
braking in the VLM regime occ urs on timescale s of a 
'few hundred Myrs'. Consistently, iTerndrup et al.1 (|200Ch 
find a value of r = 246 ± 55 Myr from a comparison of 
vsini in Pleiades and Hyades. Several studies point to 
spin-down timescales in the range of or larger than 1 Gyr 
l|Delfosse et al.l Il998l : ISills et alj I2OO0I : iBarnesI 120031 ). All 
these estimates are in order-of-magnitude agreement with 
our new constraint. Taken together, the available rotational 
data for VLM objects favours the occurence of weak rota- 
tional braking on the main-sequence, with exponential spin- 

3 Based on the currently available data, it appears that the pe- 
riod ranges in M34 and Praesepe are indistinguishable. If this 
is indeed the case, most of the rotational braking in the VLM 
range would occur at ages < 200 Myr. The spindown would be 
more rapid until about the age of M34 and would slow down con- 
siderably after that. While such a time-dependent spindown is 
certainly an interesting scenario, we do not put too much trust 
in this interpretation, because low number statistics and period 
incompleteness in Praesepe might easily be responsible for the 
effect. In addition, the mass ranges for the objects with periods 
differ slightly from cluster to cluster, which in combination with 
the strong period-mass relation (see Sect. 14. It can lead to biases. 
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down timescales larger than 200 Myr, maybe as long as a 

few Gyrs. 

According to iBarne the spin-down timescales 

for fast rotating G-stars in young open clusters are ~ 
30 Myr. Thus, among VLM objects, the spin-down happens 
at a much slower rate than for solar-mass stars: the braking 
timescales are 10-100 times longer. This results in fast ro- 
tating main-sequence objects, as already evidenced by the 
rotation vs. mass plot in Sect. 14. H Both findings, the steep 
drop in rotation periods and the similarly steep increase in 
the spin-down timescales, are complementary observational 
manifestation of a fundamental change in the angular mo- 
mentum regulation in the VLM regime. 



5 SUMMARY AND OUTLOOK 

Rotation periods have been measured from photometric 
monitoring for five stars in the open cluster Praesepe (age 
~ 750 Myr), all with masses < 0.5 Mq. Our work demon- 
strates that it is possible to obtain reliable periods for faint 
objects at the bottom of the main-sequence using wide-field 
imagers at medium-sized telescopes, as long as enough ob- 
serving time is available to provide a large number of dat- 
apoints and thus a high level of redundancy. We show that 
main-sequence open clusters like Praesepe are ideal environ- 
ments to probe the long-term angular momentum evolution 
and the underlying regulation mechanisms. 

The five periods range from 5h to almost 3.5 days and 
have been confirmed by various independent period search 
procedures as well as plausibility checks. We attribute these 
photometric variations to cool, magnetically induced spots, 
co-rotating with the objects - thus the periods likely corre- 
spond to the rotation periods. Comparing the small sample 
of periods in Praesepe with v sin i data for the coeval Hyades 
and with periods measured in younger clusters, we conclude 
that they give a reasonable first glance on the very low mass 
(VLM) period distribution in Praesepe. Particularly, it is 
unlikely that Praesepe harbours a large, undetected popula- 
tion of VLM objects with periods significantly longer than 
4 days. Still, the big caveat in the discussion of the periods 
is the small number of datapoints. As long as we do not 
have a more substantial dataset, we refrain from carrying 
out a vigorous statistical analysis and rely instead on more 
qualitative assessments. 

We find a trend of decreasing period with decreasing 
mass in the VLM regime, confirming previous claims in 
younger clusters. This trend is roughly linear in our small 
sample with P m 240 (M/Mq) h. To probe the period-mass 
relation over a broad range of stellar masses, we combine our 
dataset with periods in the coeval Hyades from the litera- 
ture. We find a dramatic change in the periods at spectral 
type K8-M2, corresponding to masses of 0.4-0.6 Mq: With- 
out clear transition regime, the periods drop from ~ 10 days 
in mid K-stars to < 4 days in early M-stars. Even consid- 
ering that the VLM period sample may be incomplete, the 
mere existence of the fast rotating M-dwarfs in Praesepe in- 
dicates a significant change in the rotational regulation. We 
note that this change coincides with the onset of chromo- 
spheric activity occuring at spectral type Ml in Praesepe 
jKafka fc HonevcutdlioOQ ) . 

Comparing the periods in Praesepe with rotational data 



in younger clusters (ages of 100-200 Myr), we find that some 
angular momentum loss is likely to occur in the VLM range. 
The exponential timescale of spin-down due to stellar winds 
is clearly longer than 200 Myr, most likely between 400 and 
800 Myr, and thus 10-100 times longer than in solar-mass 
stars. Hence, rotational braking on the main-sequence is 
clearly less efficient in the VLM regime, resulting in fast 
rotators. 

Thus, the periods in Praesepe provide evidence for two 
observational manifestations of a fundamental change in the 
rotational regulation: a clear drop of the rotation periods at 
0.4-0.6 Mq and very long spin-down timescales i n the VLM 
regim e. This is in line with previous findings : mm et al.l 
(2000) made an attempt to model the rotational evolution 
of VLM objects and conclude that it is impossible to 'simul- 
taneously reproduce the observed stellar spin-down in the 
0.6 — 1.1 Mq range and for stars between 0.1 and 0.6 Mq', 
implying that the solar-type rotational models are not ap- 
plicable at very low mass es. In the evolutionary scheme 
proposed by IBarne s (2003), the break between 'I- and C- 
sequence' (which is identical to the break in the period-mass 
relation reported in this work) can be seen as another way to 
describe a transition between two different rotational regu- 
lation regimes. In summary, the solar-type rotational evolu- 
tion picture of a transition dynamo driving a magnetic wind, 
which leads to Skumanich-type angular momentum losses, 
does not appear to be applicable in the VLM regime. 

The underlying physical reason for this fundamental 
breakdown of the solar-type rotational braking is not clear. 
We argue that it might be possible to explain it with the 
disappearance of the radiative core and the resulting inabil- 
ity to operate a solar-type interface dynamo. However, this 
transition is thought to occur at masses 0.3-0.4 Mq in solar- 
metallicity stars, i.e at somewhat lower masses than the drop 
in rotation periods. The metallicity difference between the 
Sun and young clusters like Praesepe might explain parts of 
this discrepancy. Fundamental changes in the magnetic field 
structure, the wind properties, or dynamo action (possibly 
independent of interior structure) have to be considered as 
alternative ways to understand the observational findings. 

The obvious signature of a fundamental change seen in 
the rotation periods holds great potential to use rotation 
as a probe for magnetic properties and/or interior struc- 
ture. Follow-up studies aimed a) to enlarge the rotational 
database for VLM objects on the main-sequence and b) to 
provide complementary activity-related data in this mass 
regime are thus highly desirable to clarify the open ques- 
tions. 
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